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Abstract: We present a QCD analysis of the cumulative out-of-event-plane momen- 
tum distribution in DIS process with emission of high p t jets. We derive the all-order 
resummed result to next-to-leading accuracy and estimate the leading power correc- 
tion. We aim at the same level of accuracy which, in e + e~ annihilation, seems to be 
sufficient for making predictions. As is typical of multi-jet observables, the distribu- 
tion depends on the geometry of the event and the underlying colour structure. This 
result should provide a powerful method to study QCD dynamics, in particular to 
constrain the parton distribution functions, to measure the running coupling and to 
search for genuine non-perturbative effects. 
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1. Introduction 



The success of the QCD description of e + e~ event-shape variables makes these ob- 
servables useful tools to study features of QCD radiation [1], to measure the running 
coupling [2] and to search for genuine non-perturbative effects [3, 4, 5]. 

The standard QCD analysis of event shapes involves resummations of all pertur- 
bative (PT) terms which are double (DL) and single (SL) logarithmically enhanced 
and matching with exact fixed order PT results. In addition, to make quantitative 
predictions, one needs to consider also non-perturbative (NP) power corrections. 
These standards have been already achieved for a number of 2-jet event shapes 
(1-T,C,M 2 /Q 2 ,B in e+e" annihilation [1, 6, 7] and 1-T,B in DIS [8, 9]), i.e. for 
observables which vanish in the limit of two narrow jets. 

Only recently has the attention moved to three-jet shapes in e + e~ annihilation 
(thrust minor T m [10] and D-parameter [11]). These three-jet observables exhibit a 
rich geometry dependent structure due to the fact that they are sensitive to large 
angle soft emission (intra-jet radiation). These results have been extended to jet 
production in hadron collisions [12]. The main difference between processes with or 
without hadrons in the initial state is that in the first case jet-shape distributions are 
not collinear and infrared safe (CIS) quantities, but are finite only after factorizing 
collinear singular contributions from initial state radiation (giving rise to incoming 
parton distributions at the appropriate hard scale). 

In this paper we consider the DIS proton-electron process 

e + p — > e + jets + . . . (1.1) 

in which we select events with high p t jets with p t ~ Q. The dots represent the 
initial state jet and intra-jet hadrons. This process involves (at least) three jets: two 
large p t jets (generated by two hard partons in the final state recoiling one against 
the other) and the initial state jet (generated by the incoming parton). For p t ~ Q, 
the exchanged boson of momentum q, with Q 2 = —q 2 , can be treated as elementary. 
The observable we study is 

K ut = ^\ P r\- (i.2) 

h 

To avoid measurements in the beam region, the sum indicated by Y^'h extends over 
all hadrons not in the beam direction. Here p£ ut is the out-of-plane momentum of the 
hadron h with the event plane defined as the plane formed by the proton momentum 
P in the Breit frame and the unit vector n which enters the definition of thrust major 

1 \ -v _ - 

T M = max — y \p h ■ n\ , n-P — 0. (1.3) 
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The jet events we want to select with p t ~ Q have Tm = 0(1). To select these 
events we prefer to use, instead of Tm, the (2 + l)-jet resolution variable 1/2 defined 
by the k t jet clustering algorithm [13] (see later). (By (2 + l)-jet we mean two 
hard outgoing jets in addition to the beam jet.) To avoid small and large T M we 
require y_ < 1/2 < y+ with y± some fixed values. Using this variable we have fewer 
hadronization corrections, see [4, 11]. 

The observable K out is similar to the out-of-plane jet shape studied in e + e~ anni- 
hilation [10] and in Zq production in hadron collisions [12]. In the first case the event 
plane is defined by the thrust and the thrust major axes. In the second case the event 
plane is fixed by the beam axis and the Z momentum. The reason to analyse distri- 
butions in the out-of-plane momentum is that the observable is sufficiently inclusive 
to allow an analytical study. The analysis of the in-plane momentum components is 
more involved since one needs to start by fixing the jet rapidities. 

Our analysis of K out will make use of the methods introduced for the study of 
the two observables in [10] and [12]. In the present case we have one hadron in the 
initial state, so that incoming radiation contributes both to the observable K out and 
to the parton density evolution. To factorize these two contributions we follow the 
method used in hadron-hadron collisions [12]. The result is expressed in terms of the 
following factorized pieces: 

• incoming parton densities obtained by resumming all terms a™ ln n fi/K out (// is 
the small factorization scale needed to subtract the collinear singularities and 
K out is the hard scale for this distribution); 

• "radiation factor" characteristic of our observable. Its logarithm is obtained by 
resumming all DL and SL terms (a™ ln™ +1 K out /Q and a™ ln n K out /Q respec- 
tively). 

• matching of the resummed result with the fixed order exact calculations [14, 15]. 

The factorization of the incoming parton densities and the radiation factor is the 
crucial step for the present analysis. This result is due to coherence and real/virtual 
cancellations (see [12] and later). As a result, after this factorization procedure, the 
radiation factor is a CIS quantity similar to the ones entering in the e + e~ observables. 

In order to make quantitative predictions we need to add to the above PT result 
the 1/Q-power corrections. We follow the procedure successfully used in the analysis 
of e + e~ jet-shape distributions [4]. The definition of the event plane makes our 
observable sensitive to hard parton recoil. Here only the two final state hard partons 
can take a recoil, while the initial state one is fixed along the beam axis. This 
simplifies the treatment both of the PT distribution and of the interplay between 
PT and NP effects, which is a characteristic feature of all rapidity independent 
observables. 
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The 1/Q coefficient is expressed in terms of a single parameter, cto (/•*/), given 
by the integral of the QCD coupling over the region of small momenta kt < Hi 
(the infrared scale /// is conventionally chosen to be /ij — 2 GeV, but the results 
are independent of its specific value). Effects of the non-inclusiveness of K out are 
included by taking into account the Milan factor M. introduced in [6] and analytically 
computed in [16]. The NP parameter a (fij) is expected to be the same for all jet 
shape observables linear in the transverse momentum of the emitted hadrons. It 
has been measured only for 2-jet event shapes and appears to be universal with a 
reasonable accuracy [5]. It is interesting to investigate if this universality pattern 
also holds for near-to-planar 3-jet observables. 

In order to improve the readability of the paper, in the main text we present 
only the results and discuss their physical meaning, leaving the detailed derivation 
of the results to a few technical appendices. In section 2 we define the distribution 
and specify the phase space region of K out in which we perform the QCD study. In 
section 3 we describe the PT and NP result for the K out distribution. We stress how 
the answer has a transparent interpretation based on simple QCD (and kinematical) 
considerations. In section 4 we improve our theoretical prediction by performing the 
first-order matching and present some numerical results. Finally, section 5 contains 
a summary, discussion and conclusions. 



2. The process and the observable 

We work in the Breit frame 

g=|(0,0,0,2), P= ^-(1,0, 0,-1), *B = ^y, (2.1) 

in which 2x B P + q is at rest. Here P and q are the momenta of the incoming proton 
and the exchanged vector boson (7 or Z ). 

In this frame, the rapidity (with respect to the direction of the incoming proton) 
of an emitted hadron with momentum ph is given by 

l ln ( 1+ « ). (2. 2) 

To avoid measurements in the beam region, the outgoing hadrons ph are taken in the 
rapidity range 

Vh < Vo - -Infancy , (2.3) 

which corresponds to a cut of angle O around the beam direction. Similarly, the 
sums Yl'h m (1-2) an d (1-3) extend over all hadrons with rapidities in the range (2.3). 
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To select jet events with p t ~ Q we use the (2 + l)-jet resolution variable y 2 
introduced in the /c r algorithm for DIS processes [13]. For completeness we recall its 
definition. 

Given the set of all outgoing momenta one defines the "distance" of ph from the 
incoming proton momentum 

VHP = -^E 2 h (l- cos 9 hP ). (2.4) 

For any pair p^ and p^' of outgoing momenta one also defines the "distance" 

2 

Vh'h" = q5 min{££,,^,,}(l - cos9 h , h „) . (2.5) 

If Dhp is smaller than all i/h'h", the hadron p^ is considered part of the beam jet 
and removed from the outgoing momentum set. Otherwise, the pair of momenta 
Ph',Ph" with the minimum distance ywh" are substituted with the pseudo-particle 
(jet) momentum p h m = p h / + Ph" (-E-scheme). The procedure is repeated with the 
new momentum set until only two outgoing momenta Phi,Ph 2 are left. Then the final 
value of y 2 is defined as 

y 2 = mm{y hlP ,y h2 p,y hlh2 } . (2.6) 
To select jet events with p t ~ Q, as stated before, we require 

y_<y 2 <y + , (2.7) 

with y± fixed limits. 

The distribution we study is then defined as 

da(y±,K out ) ^ f da, 



dx-s dQ 2 



£ / iBq^-^-^ q k.-X>ri . (2.8) 

m " ^ \ h=l / 



with da rn / dx^dQ 2 the distribution for m emitted hadrons in the process under con- 
sideration. Considering the cross section for the y 2 -range (2.7) 

^-E/^e fo+ - 92 )e te -„_), ( , 9) 

we have the normalized distribution 

Due to the rapidity limitation (2.3) in the definition of y 2 , K out and the event plane, 
the distribution will depend on i] . To avoid a strong ^-dependence we will consider 
i]o and K out in the range (see later) 

^out > Qe-». (2.11) 
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The fact that this distribution is rich in information on the hard process is clear from 
the fact that it depends not only on the observable K out but also on the variables 
y± which define the geometry of the jet events with p t ~ Q. In the range (2.11) the 
PT result does not depend on i] , to our accuracy. However, as we shall discuss, the 
power correction depends linearly on i] . This is due to the fact that the contributions 
to the observable are uniform in rapidity. 



3. QCD result 

The QCD calculation of the distribution (2.10) is based on the factorization [17] of 
parton processes into the following subprocesses: 

• elementary hard distribution; 

• incoming parton distribution; 

• radiation factor corresponding to the observable K ont . 
These factors are described in the following subsections. 

3.1 Elementary hard process 

For the elementary hard vertex 

qP 1 ^P 2 P 3 , (3.1) 
we introduce the kinematical variables (see Appendix A) 

£- (PlP2) x- Q2 >x n (3 2) 

The invariant masses are 

Ql b = 2(P a P h ), Q\ 2 = i -Q\ Q\,= l -^Q\ Ql,= l -^Q 2 - (3.3) 

The thrust major for the process (3.1) is given by Tm = 2^£(l — £)(l — x)/x. The 
substitution £ — > (1 — £) interchanges P 2 and P 3 , so we distinguish P 2 from P 3 by 
assuming 

(P X P 2 ) < (P X P 3 ) , (3.4) 

which restricts us to the region < £ < \. The variable y 2 for the elementary vertex 
is given in terms of the variables in (3.2) by 

. 2 , ^ e \+£ 



x ' 1 + 2 ^' ^ 

- \ /l _ ^ K . ( i £)(l , r) for ^ > 1 + e ' 



x 



x{l-t) + t{l-x) ' l + 2£ 
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Inverting this one finds £ = ^(x,y 2 ). This function and the relative phase space are 
discussed in Appendix A. 

We consider now the nature of the involved hard primary partons. We iden- 
tify the incoming parton of momentum Pi by the index r = q,q,g. Since (3.4) 
distinguishes P 2 from P 3 , in order to completely fix the configurations of the three 
primary partons, we need to give an additional index 5 — 1, 2, 3 identifying the gluon. 
Therefore the primary partons with momenta {P 1 , P 2 , P 3 } are in the following five 
configurations 

r = g, 5 = 1 -> {gqq} , or {gqq} 
r = q, 5 = 2 -> {qgq} , 

r = q,5 = 2 - {qgq}, (3-6) 
r = q, 5 = 3 -> {gg#} , 
r = g, 5 = 3 -> {ggc/} . 

In Appendix A we give the corresponding five elementary distributions da Tt sj, with 
/ the fermion flavours. The presence of the index r as well as 5 is due to the 
parity-violating term in the cross-section associated with Z Q exchange, such that 
the elementary cross-sections differ for incoming quark and antiquark of the same 
flavour. If we consider only photon exchange then the index r is redundant. 

3.2 Factorized QCD distribution 

The process (1.1) is described in QCD by one incoming parton of momentum p\ (in- 
side the proton), and two outgoing hard partons p 2 ,P3 accompanied by an ensemble 
of secondary partons h L 

qpi -> P2P3h ■ ■ -k n . (3.7) 

Taking a small subtraction scale [i (smaller than any other scale in the problem), we 
assume that p 1 (and the spectators) are parallel to the incoming hadron, 

Pl = Xl P. (3.8) 

Therefore, the observable we study is 

^out = \P2x\ + \P3x\ + ^2\hx\, (3.9) 

i 

where the x-axis corresponds to the out-of-plane direction, the z-axis to the Breit 
direction, and the event plane is the y-z plane. For large y 2 the hard partons p 2 ,p3 
are emitted at large angle. For small i^out) P2,P3 and the secondary parton momenta 
ki are near the event plane. 
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The event plane definition (1.3) corresponds to the condition 



P2x + P3x +J2 ki *' ( 3 - 10 ) 

U D 

which, together with momentum conservation, leads directly to 

P2x = ~ / J ki x — 2 / ki x , P3x = ~ / ^ix ~~ 2 / > ' (3-11) 

Here, by [/ and .D we indicate the up- and down-regions corresponding to partons hi 
with kiy > and A;^ < respectively. Again, by ^2 we indicate the sum restricted 
to secondary partons in the region (2.3). By Yl" we indicate the sum restricted to 
secondary partons ki in the "beam- region" with r\i > i] . 

We perform the QCD analysis at the accuracy required to make a quantitative 
prediction: DL and SL resummation, matching with exact fixed order results, and 
leading 1/Q power correction. The analysis is similar to the one performed in [10, 
11, 12]: the starting point is the elementary process (3.1). Then one considers the 
secondary radiation (soft and/or collinear) in a 3-jet environment. Finally one takes 
into account the exact matrix element corrections and power corrections. 

The application to the present case is described in detail in Appendix B where 
we show that the distribution can be expressed in the following factorized structure 

da(K out ,y ± ) ^ f XM dx /*+ / da T 



dx B dQ 2 TiS j Jx B 



(3.12) 



where xm is a function of y_, the limits £± = £±(x, y±) select jet events with p t ~ Q 
and da are the distributions for the elementary hard process (3.1). See Appendix A 
for the elementary distributions and kinematics. 

The distributions X, which resum higher order QCD emission, can be expressed 
in the following factorized form 

Zr,8j = C TtS {a s ) ■ V Tj (—, K ont j ■ A 5 {K out , x, f , Q, i] ) . (3.13) 

Here we describe the various factors: 

• the factor V T j is the incoming parton distribution. It is given, for the various 
cases, by the quark, antiquark or gluon distribution inside the proton 

V qJ (x, K out )=q f (x, K out ) , V qJ (x, K out )=q f (x, K out ) , 
V gJ (x, K out ) =g(x, K out ) . 

We show that here the hard scale is fixed at K out and that the dependence 
on r/o can be neglected as long as we take K out sufficiently large in the range 
(2.11); 
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• the distribution As is the CIS radiation factor which resums powers of In K out / Q 
and is a CIS quantity. It is sensitive only to QCD radiation and therefore does 
not depend on the flavour (we neglect quark masses). There are various hard 
scales in As (given in terms of the Q 2 ab in (3.3)) which are determined by the 
SL accuracy analysis. This quantity is similar to 3-jet shape distributions one 
encounters in e + e~ annihilation processes [10, 11] and in hadron collisions [12]; 

• the first factor is the non-logarithmic coefficient function with the expansion 

C T>s (a s ) = l + Cl ^ + (a 2 s ) , a s = a s {Q) . (3.15) 

It takes into account hard corrections not included in the other two factors and 
is obtained from the exact fixed order results. 

The factorization of the first two pieces is based on the fact that contributions to 
A (to V) come from radiation at angles larger (smaller) than K out /Q. This implies 
that one is able to reconstruct the parton densities V as for the DIS total cross 
sections in which one does not analyse the emitted radiation. The only difference 
here is that the parton density hard scale is given by -Kout while in the fully inclusive 
case of DIS total cross section the hard scale is Q. As a result of this factorization, 
the radiation factor A can be analysed by the same methods used in e + e~. 

The basis for the factorized result (3.13) has been discussed in [12] and it is 
based on coherence and real/virtual cancellations. Since this is a crucial point for 
our analysis we recall in some detail the relevant steps. 

• Since K out is a CIS (global [18]) quantity (its value does not change if one of 
the emitted particles branches into collinear particles or undergoes soft brems- 
strahlung), within SL accuracy we can systematically integrate over the final 
state collinear branchings. 

• Factorizing the phase space and the observable (see (B.4)) by Mellin and Fourier 
transforms one has that each secondary parton of momentum ki contributes 
with an inclusive factor 

1(h) = 1 - e(h) ■ U{h) , (3.16) 

where 1 is the virtual term while the factors e and U are the contribution from 
the real emission. The factor e(ki) accounts for the energy loss of the primary 
incoming parton pi due to emission of collinear secondary partons. Therefore, 
for ki non-collinear to p\ we have e(fcj) = 1, while for ki collinear to p\ we have 
e(fcj) = zf~ l with 1 — Zi the energy fraction of ki with respect to p\ and iV the 
usual Mellin moment of the anomalous dimension, see (B.7). The factor U(ki) 
depends on the observable and in our case is given by (B.8). 
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• The crucial point which is the basis of the factorization in (3.13) is that, to SL 
accuracy for K ont <C Q, we can replace (see (B.16)) 

[1 - U{h)} ~ Q(k u - ho) , k ~ K ont , (3.17) 

so that we can write 

I(k t ) ~ (hi-ko) + (*o - • (3.18) 

The first term contributes to the radiation factor while the second reconstructs 
the anomalous dimensions for the various channels. The two contributions do 
not interfere and, as discussed above, one obtains the factorized expression 
(3.13). 

In the next two subsections we will describe our result for the radiation factor As- 
First we describe the PT contribution obtained by resumming the logarithmically 
enhanced terms at SL accuracy. Then we describe the leading NP corrections orig- 
inating from the fact that the (virtual) momentum in the argument of the running 
coupling cannot be prevented from vanishing, even in hard distributions. Matching 
with the exact first order result is considered later. 

3.3 The PT radiation factor and distribution 

The PT radiation factor A^ T can be expressed as (see Appendix B) 

A PT = e -Rs(K-^,Q) . S f as ln (3 ig) 

V ^out / 

It resums DL and SL contributions originating from the emission of soft or collinear 
secondary partons hi in (3.7). The PT radiator R$ in the first factor contains the DL 
resummation together with SL contributions coming from the running coupling and 
the proper hard scales. The second factor is the SL function which resums effects due 
to multiple radiation, including hard parton recoil. For K out in the region (2.11), this 
PT result does not depend on i] . The analysis leading to this result is very similar 
to that for other 3-jet distributions [10, 11, 12] (see Appendix C). Here we report 
the relevant expressions and illustrate their physical aspects. 

The PT radiator R$ is given by a sum of three contributions associated with the 
emission from the three primary partons P a . For the configuration 5, it is given by 

Rs = E^MiCUW) , (3-20) 

a=l 

where 

Qa ^ = Qb = Qab i 

cf = cf = C F , 



QasQ 



5b 



(S) 



.(8) 



Qab 

Ca , 

go 

P 4JV C 



(3-21) 
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with indices a, b corresponding to the fermions and 5 to the gluon, i.e. a ^ b ^ 5. 

As is typical for a 3-jet quantity, the scale for the quark or antiquark terms 
of the radiator is the fermion invariant mass, while the scale for the gluon term is 
given by the gluon transverse momentum with respect to the fermion system. The 
rescaling factors ffi take into account SL contributions from non-soft secondary 
partons collinear to the primary partons P a . 

In (3.20) r is the following DL function 

r(K2,C?) S f f^.n|, (3,2) 

and k is the out-of-plane component of momentum and a s is taken in the physical 
scheme [19]. The rescaling factor 2 in the running coupling comes from the integration 
over the in-plane momentum component. The fact that K out is the lower bound in 
the PT radiator, to this order, comes from real/virtual cancellation (see (3.17) and 
(B.16)). The exact expression of the hard scales and the rescaling factors in (3.21) 
is relevant at SL level. 

The factor S is expressed in terms of the SL function 

r' = r'(K-^ Q) = ln ® , (3 .23) 

A out 

given by the logarithmic derivative of r, apart from terms beyond SL accuracy. We 
find 



g-7B Ct r' 



r(i + C T r f ) 



- ■ T (cS r') T (cg> r') , C T = 2C F + C A , (3.24) 



where Ct is the total colour charge of the primary partons (see (3.1)). The first factor 
is the same for all 3-jet shape emission distributions. The second factor depends on 
the specific observable. The function T is given by 

r(l±2) 

It takes into account the correct kinematics for the emission and in particular the 
effect of the recoil of the two outgoing primary partons in (3.1) due to the emission 
of secondary QCD radiation. Here C[ S J is given by the colour charge of the out- 
going parton a = 2, 3 plus half of the charge of the incoming parton, in the given 
configuration, 

Cg = \C? + Cf . (3.26) 

This combination of colour charges enters due to the kinematics that defines the 
event plane (see (3.11)), which leads to the vanishing of the sum of the out-of- 
plane momenta in the "up-region" , i.e. with positive components, and also in the 
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"down-region" . (We are working in the regime in which the effect of the rapidity cut 
is negligible.) The momenta of outgoing primary partons #2 and #3 are in the up- 
and down-region respectively. The incoming parton $4 is instead along the Breit 
axis and emits equally into both regions. Notice that Cf^ + = C T . 

Finally the DL and SL resummed PT part of the distribution J is obtained from 
(3.13) by using the radiation factor in (3.19) 

2*J, = C T , s (a s ) ■ V Tj (^, Kj) ■ Ar (K out , x, £, Q) . (3.27) 

The coefficient function will be computed at one-loop by using the numerical program 
DISENT of Ref. [14] and subtraction of the one-loop contribution already contained 
in the two factors V and A. 

To first order in the PT expansion we have 

d\nK out TT f-f K out K sJ 

a=l 

The absence of the factor \ in the logarithm compared with that of equation (3.22) 
is due to the event-plane kinematics (see (3.11)). For a single gluon k x (emitted 
from pi,p2 or p 3 ) only one of the two hard partons P2 or p 3 takes recoil equal to 

\Pax\ = \k x \ = K on t/2. 

3.4 The distribution including NP corrections 

As in other cases of jet-shape distributions, the leading NP correction corresponds to 
a shift in the PT distribution. This is due to the fact that, in the Mellin representation 
of X, the PT radiator is affected by a NP correction with leading term linear in the 
Mellin variable. Thus its effect corresponds to a shift in the conjugate variable, K out , 
and one has 

Zr A /(*c«t) = J 3/ (*o«t - • (3-29) 

The quantity 8K^ t , which can be cast in the following form 

f r {S) 3 r (5) 3 i 

SK& = A NP I C -g ■ E ^ + ■ E y« I , (3.30) 

corresponds to the integral over the infrared region of the soft gluon distribution 
with weight \ky\, the contribution to K ont from the very soft gluon generating the NP 
contribution. The first factor A NP corresponds to the momentum integral (including 
\k x \ and the running coupling), while is the rapidity interval (more precisely 
the logarithmic integral of the angle which the very soft gluon forms with the event 
plane). 
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The parameter A NP , given in (D.6), is expressed in terms of the NP parameter 
ao{fii) and the Milan factor M. [6, 16] which takes into account effects of the non- 
inclusiveness of K out . A NP also contains renormalon cancellation terms. The quantity 
a (/ij) is the integral of the running coupling in the infrared region 



dkt 



<xoM = I —a s (k t ) . (3.31) 
Jo Hi 

The infrared scale /// is conventionally chosen to be /// = 2 GeV, but the results are 
independent of its specific value. Both ao(fij) and M. are the same for all jet shape 
observables linear in the transverse momentum of emitted hadrons. 

The derivation of the rapidity intervals is reported in Appendix B, see 
(B.39). The expressions for 6=1 are simply given by 

y 2 '?=^»=ln^. (3.32) 
For the diagonal cases a = b = 2, 3 we have the behaviour 

F 2 ? - ^ ^ + O (r') , F3? - ^ ^ + O (r') , (3.33) 

-ft- out out 

with r' given in (3.23), Q*)p the hard scales in (3.21) and ( = 2e~ 2 the rescaling 
factor which accounts for the mismatch in the integration over an angle between two 
vectors and a vector and a plane. The expressions for the off-diagonal cases have more 
complex behaviours. We have to distinguish two regions. For a s In 2 Q/K out < 1 we 
have (see (B.42)) 



1 23 — / ) J 32 — 



I ) *32 — ; ) 

with corrections of 0(1). For a s In 2 Q/K out 3> 1, we have (see (B.41)) 



(3.34) 



Y 3 f * In ^ + O (r') , c In ^ + O (r') . (3.35) 

-"■out ^ out 

Before illustrating and discussing these expressions, we observe the following 
features of the distribution obtained in (3.29). Since SK^ t depends on K out , the 
PT distribution is actually deformed. The dependence on 770, the available rapidity 
interval for the measurement, enters only through the NP shift (the PT distribution 
is independent of 770, to our accuracy, as long as we take i^out m the region (2.11)). 
The dependence of ^i^out on 2/2 (through Q^), implies that the deformation of the 
PT distribution differs for different geometries of the process. To this order (leading 
1/Q power) we can neglect NP corrections arising from the anomalous dimension 
[4] which are of subleading order 1/Q 2 . Finally, we observe that the presence of 
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lnK out and contributions to the shift are features common to the case of 

rapidity independent observables such as, broadening [7], T m [10] or the out-of-plane 
momentum in hadronic collisions studied in [12]. 

We now discuss how the above expressions for the rapidity intervals are 
based on simple QCD considerations and on the kinematics of the event-plane (see 
(3.11)). In the following we focus on the term (3.30) 

^ • ■ Yg , (3.36) 

the other is obtain by exchanging 2^3. It corresponds to the contribution in which 
a PT gluon is emitted in the up-region (this happens with probability C[ 2 /Ct) and 
a very soft gluon, generating the NP correction, is emitted off the primary parton 
6 = 1,2,3 (with colour factor Cfp). For this contribution we consider the various 
cases. 

Case 6=1. The incoming parton pi, emitting the very soft gluon, is never af- 
fected by recoil since it is aligned along the Breit axis. Therefore the rapidity interval 
Y" 2 i is fixed by the boundary r] , the available rapidity interval for the measurement 
(see (2.3)). The remaining term InxQ^/Q is related to the standard parton #1 
hard scale and a boost from the Breit frame of full process to the hard elementary 
vertex, see (D.10); 

Case 6 = 2. The outgoing parton p 2 and the PT gluon are in the same region 
(the up-region) so that, from the event plane kinematics (3.11), P2 undergoes recoil 
with \p 2x \ = K out /2. Since \n((Q 2 S ^ /\p2x\) fi xes the rapidity interval for the very soft 
gluon emitted by p 2 , we obtain the result (3.33); 

Case 6 = 3. Here we cannot be limited to considering the contribution from 
a single PT gluon emitted in the up-region since, from the event plane kinematics 
(3.11), one would have p 3x = and the rapidity interval Y 2 3 would diverge (for a 
zero mass gluon). We then need to consider high-order PT contributions which push 
parton #3 off the plane by a certain p 3x . The distribution in the recoil momentum 
P3 X has two different regimes: 

• for a s In 2 Q/K out <C 1 there are few secondary partons and the PT distribution 
is given by Sudakov form factors. Integrating ln|p3 X |, the rapidity interval 
contribution, we obtain the behaviour in (3.34) (here is the coefficient of 
the DL exponent of the \p$ x \ Sudakov factor); 

• for a s In 2 Q / K out 1 the PT radiation is well developed, so that pz takes recoil 
with \pz x \ ~ -^out and we obtain the behaviour (3.35). 

4. Matching and numerical analysis 

The final step needed to obtain a quantitative prediction is the calculation of the 
non- logarithmic coefficient function (3.15) from the exact matrix element results. We 
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use the numerical program DISENT of Ref. [14] , which includes only the 7 exchanged 
contribution, and we compute the first order coefficient c\ by subtracting the DL and 
SL terms already contained in the PT resummed result. The Z contribution will be 
important at large values of Q. For example, for Q = 20 GeV the Z contribution to 
the cross section is of order 0.5% (see (A. 8)). 

Taking into account only 7 exchange simplifies considerably the expressions since 
we only need 5 to identify the configuration, while r becomes redundant. The dis- 
tribution (2.8) and the cross-section for the high p t jet (2.9) are given by 

dx B dQ 2 jjJ XB x 4 \dxdtdQ* )''^\x^ 

(4.1) 

where xm and £± = £±(x,y±) are found in Appendix A. The hard elementary 
distribution is given by (see Appendix A) 

(Jtyil) 2 

x 2 +e 



;i-^)(i-o 



+ 2(l+x0 



C£ = C F • Axi , Ci = 8z(l-z) , <%, L {x,Z)=C* /L (x,l-Z) , 

where y = Q 2 /(x B s) with ^/s the centre of mass energy of the process. 

We start with the resummed PT result X PT (3.27) which in this case takes the 
form 

with the incoming parton distribution given by 

V$(x,K mt ) = \ 9 \' ° Ut) (4.4) 

The PT radiation factor A^ T is the same whether we include the Z or not, and is 
given in (3.19). 



4.1 Matching 

the first or 

of the coefficient function (see (3.15)), we expand the normalized resummed 



We now discuss the first order matching. To determine the first order term 
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distribution to one loop 



dx B dQ 2 { dx B dQ 2 j ^ ^ 

= 1 + ^ (ft* In 2 ^- + G<? (!/± ) In -5- + ^ {y± , K M )) + ■■■ 

27T \ A out A out / 

and compare with the result provided by the numerical program DISENT [14]. From 
the resummed result we have G i2 = —2Ct and 



5,/ 



(4-6) 

4 ^ 6a ln Q a s (g)91ngj ^U' g J • 

We have verified that the above coefficients G12 and G$ are correctly reproduced 
by DISENT and by subtraction we deduce the non-logarithmic one loop coefficient 
c?(y ± ,K out ). 

We use the first-order (Log R)-matching prescription and write 



mat 



e gci(y±,K out ) ^ ^ ^ 



with a s = «Ms(g)- We also implement the correct kinematical bound for K out < K^ t 
(with K™ t obtained from DISENT) by replacing K out by K' out , implicitly defined by 

-f^out xl out zv out -f^out 

in the final PT result (4.3) so that the differential distribution vanishes for -Kout = 
K^ t . Finally we make the replacement (3.29), (3.30) in order to include the NP 
correction. 



4.2 Numerical results 

To summarize, the QCD prediction for the normalized distribution £(y±,if ut) in 
(2.10) is obtained (for exchanged 7) from the ratio of the two distributions in (4.1) 
in which the integrand X is given by (4.3) with the coefficient function given in (4.7) 
and using the substitution (4.8) of K out . The radiation factor A is given, at PT level, 
by (3.19). Leading power corrections are included by performing the substitutions 
in (3.29) and (3.30). 

We now report some plots for cE(y±, K out )/dK out (data are not yet available) 
for s = 98400 GeV 2 and the rapidity cut r] = 3. A smaller rapidity cut, i.e. a larger 
value of 770, would have two advantages: the PT resummation is valid for smaller 
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Figure 1: QCD prediction for x B = 0.1, Q 2 = 400 GeV 2 and y_ = 0.1. In the NP 
correction we have taken ao(2GeV) = 0.52. 
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Figure 2: The same as Fig.l but with Q 2 = 900 GeV 2 . 



K out , and the NP shift becomes larger; but this needs to be balanced against the 
experimental resolution and the need to exclude the proton remnant. 

The QCD predictions are all obtained for the following choices: the parton den- 
sity functions in [20], set MRST2001; the NP parameter ao(yU/) = 0.52, a value in 
the range determined by the analysis of 2-jet observables in e + e~ annihilation [21] 
and the running coupling a s (Mz) = 0.119. 

Figs. 1 and 2 show the differential K ont distribution for Q = 20 GeV and Q = 
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Figure 3: The full (PT+NP) differential K on t distribution for two different values of y_, 
for fixed y + = 2.5, x B = 0.01 and Q 2 = 400 GeV 2 . 



30 GeV, and x-q — 0.1. In order to prevent Tm from being too large it is enough to let 
y + go to its maximum kinematically allowed value y + = 2.5 (see (A. 4) and (A. 5)). 
From these pictures one notices the 1/Q dependence of the power correction. One 
can also observe that the shift is larger at smaller K ont , so that the distribution is 
squeezed. This is due to the In K out dependence of the NP correction (see discussion 
in section 3.4). 

The fact that our predictions are also geometry dependent can be seen from 
Fig. 3, where we plot the PT+NP curve for two different values of ?/_, fixing y + = 2.5. 

Finally, Fig. 4 shows how the distribution depends on the value of x-q. The 
distribution dies faster by increasing x-q- Indeed, for increasing x B the centre-of- 
mass energy of the hard system decreases, so that the phase space for producing 
hadrons with large out-of-plane momentum is reduced. Moreover, by increasing 
x-q the configuration with incoming quark or antiquark (5 = 2, 3) becomes more 
probable. 



5. Discussion and conclusion 

In this paper we have extended our knowledge of three-jet physics to the case of a 
near-to-planar observable in hard electron-proton scattering, defined analogously to 
the thrust-minor in e + e~ annihilation. On the PT side, this observable exhibits a 
rich colour and geometry dependence. Large angle gluon radiation contributes in 
setting the relevant scales for the PT radiator (see (3.21)). The radiator itself is 
a sum of three contributions, one for each emitting parton, each one characterized 
by a different hard scale: for a quark or antiquark it is the invariant mass of the 
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fermion system, for a gluon it is its (invariant) transverse momentum with respect 
to the fermion pair. Such a structure, due to the universality of soft radiation, is 
found to be common to all near-to-planar 3-jet shape variables encountered so far 
(see [10, 11, 12]). 

Since the observable K ont is uniform in rapidity, hard parton recoil contributes 
both to the definition of the event plane and to the observable. It contributes then 
to the PT distribution at SL level. 

The thrust major axis is determined only by the hard parton system and cannot 
be changed (at least at SL accuracy) by secondary soft or collinear parton radiation. 
This implies that there exists no correlation between the up- and down-regions (see 
the definition of the event-plane (3.11)). This property makes the analysis of our 
observable much simpler than that of the thrust minor T m distribution in e + e~ . In 
particular, we are able to write a closed formula for the SL function S$ in (3.24) 
which embodies the contributions of hard parton recoil. 

Our observable exhibits 1/Q power corrections arising from the running of the 
coupling into the infrared region. For K out ^> Aq CD they can be taken into account as 
a shift of the PT distribution (see (3.29)). The shift 8K^} t depends logarithmically 
on the observable itself, so that the effect of the shift is actually a deformation 
of the PT distribution. Since K out accumulates contributions from partons of any 
rapidity, one has to carefully consider effective rapidity cut-off. Along the beam axis 
the rapidity is bounded by the experimental resolution rjo- Along the direction of 
the two outgoing hard partons, it is their displacement from the event plane which 
provides an effective rapidity cut-off. Averaging such a NP correction over the hard 



19 



parton PT recoil distribution gives rise to contributions of order lni^ out or 
according to the K out regimes. The form of the power corrections can be simply 
interpreted on the basis of the event plane kinematics (see (3.11)), as explained in 
detail in section 3. 

As in other 3-jet observables, the magnitude of the NP shift is expected to 
be roughly twice as large as the corresponding one for 2-jet shape variables (the 
characteristic weight 2C F of the DL contribution of the logarithm of 2-jet distribution 
becomes here a 2Cf + Ca weight). This implies that higher order power corrections 
may become important, see [22]. The comparison with data (not yet available) will 
be needed to check the validity of our method. 
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A. Kinematics and elementary partonic cross-sections 



For the elementary process (3.1) we may write the parton momenta Pi, P2 as 
P 1 = £(1,0, 0,-1), P 2 = Q(z ,0,T M ,z 3 ), 

^ (i-fl(i-s) _ (i-0(i-^) 

^0 — ? H , ^3 — ? 



(A.1) 



.7; 



in terms of the variables in (3.2). Distinguishing P 2 and P 3 according to (3.4), the 
variable £ is given in terms of T M by 



1 - \ 1 



rT 2 

1 — X 



(A.2) 



and in terms of y 2 by the inverse of (3.5), which is 



l-x 

2(2x-l) 



Ax(2x - l)|/2 



1 + 



{l-xf 



- 1 



X 2 V2 - (1 - X) 2 

(2x - 1)(1 - x + xy 2 ) ' 



3 1 n 2(1 -x) 2 

4 <x<1 ' 0<2/2< iT^T) 

3 ^ ( 2(1 - xf l-x 

4 ' x(2x — 1) ^ 2 a; 

(A.3) 



Note that at any fixed value of y 2 , the freedom to vary x results in events with a 
range of T M values. If y 2 < |, the kinematic upper limit on x is x = at which 



£ = \ and the lower bound on Tm is therefore ^fy^. The lower limit on x is the 
Bjorken variable x B , at which £ = £,{x-Q,y 2 ) < x-&y 2 . (This is obtained by expanding 
(A.3) about xb = 0: the difference from equality is of order x B .) Consequently, 
Tm < 2\/ r y~2- Therefore, restricting ourselves to events in the range (2.7), with 
y- < I < y + , leads to a selection of events whose values of T M lie in the range 



y_ < T M < 2y/y^ . 



(A.4) 

This of course does not select all events with T M in this range; it is merely a means of 
selecting events whose T M does not lie outside it, so that we do not need to consider 
logarithms of T M . The phase-space in terms of the variables (x,£) then becomes: 



Xg < X < Xm , %M 



1 + 3/- 
1 



, y_<l/3 



4y_ 



V- > 1/3 



1 

2 ' 



x > 



4y+ 



(A.5) 
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Since we consider only diagrams with the exchange of a single photon or Zq, the 
individual partonic cross-sections may be decomposed according to: 



s,f 



a a 



dxd^dQ 2 



{cf(Q 2 ) [{2-2y + y 2 )C r /(x,0 + 2(l-y)C T L 5 (x,0 
+ Df(Q*)y(2-y)C T /(x,0}, 



(A.6) 



into transverse, longitudinal and parity-violating terms. Here y is defined by 

(Pq) Q 2 



y 



(PPe 



X B S 



(A.7) 



where P e is the momentum of the incident electron, and s is the centre-of-mass 
energy squared of the collision (we neglect the proton mass) . The flavour-dependent 
functions C^(Q 2 ) and D^(Q 2 ) show how the 7 and Z exchange diagrams combine: 



C f (Q 2 ) 



2e f V f V e f Q 2 \ (Vf + A 2 )(y 2 + A 2 e ) ( Q 2 \ 



sm 2 29 w \Q 2 +M 2 / 



sin 29 w 



f 2e f A q A e ( Q 2 \ 4V f A f V e A e f 

W ' " sm 2 20 w \Q 2 + M 2 ) + sin 4 2^ \ 



Q 2 \ 



r M2J ' (a.8) 



Q 2 +M 2 J ' 



The functions C^ L ^ 3 (x,^) are the one loop QCD elementary square matrix elements 
for the hard vertex (there is no contribution of order O (a®) since we require trans- 
verse jets). We have [23] 



x 2 +e 



(l-x)(l-0 



Cf = C F -4xt, 



Cf = C f 



x 2 +e 



(l-x)(l-0 



+ 2(l + x£) 



+ 2(x + 



Cf = [x 2 + (1 - xf] 



e 



Cf' 1 = &r(l - x) , 



cf - 1 = , 



(A.9) 

B. Resumming QCD emission 

The QCD results illustrated in section 3 are based on calculations similar to those 
performed for other 3-jet shape distributions in e + e~ annihilation [10, 11] and in 
hadronic collisions [12]. Here we follow the method and use the results developed 
and obtained there. We do not report the full details of the calculations but simply 
discuss the main features characteristic of the present distribution: for details refer 
to the later appendices and to the previous papers. In particular here we 

• deduce, in the present context, the factorized structure in (3.13) and identify 
the hard scale entering the incoming parton distribution V; 
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• obtain the PT radiation factor at the SL accuracy (3.19); 

• compute the leading order NP corrections giving rise to the shift (3.29); 

• show that, for K ont in the range (2.11), the r] dependence enters only in the 
NP shift. 

B.l Resummation of the distribution 

Considering the region K out <C Q, the starting point is the factorization of the 
square amplitude for the emission of n secondary soft partons in the process (3.7) 
(contributions from hard collinear secondary emission are included later). We have 

\M T>s>f>n ( k i . . . k n )\ 2 ~ |M tA/i0 | 2 • S s , n (h ...k n ). (B.l) 

The first factor is the Born squared amplitude which gives rise to the elementary hard 
distribution da Tj sj in (3.12). The second factor is the distribution in the soft partons 
emitted from the system of the three hard partons p±,P2 and p 3 in (3.7). It depends 
on the colour charges of the emitters which are identified by the configurations 5 = 
1,2,3, the index of the primary gluon momentum. Since soft radiation is universal 
and does not change the nature of incoming parton, Ss >n does not depend on r and 
/. For soft emission, the primary partons p a differ from the hard primary parton 
momenta P a (depending on Q,x, £, see Appendix A) by small recoil components. 
By using (B.l), the soft contributions to the cross-section are resummed by 

da(K out ,y±) /*+ / da T ~- ^ ^ 



dx B dQ* ^ ./,,,, .i< 



dX ^ {dx^dQ^I l *°iPtAw) 

x£An /— fd* n -s 8 , n , 

^— ' n\ f ■» J ircji J 

1=1 



(B.2) 



with xm given in Appendix A and x\ the momentum fraction (3.8) of the incoming 
parton and \x the (small) subtraction scale. The phase space d$ n fixes the observable 
K out , the event plane and x, the Bjorken variable for the hard elementary distribution 
da . The momentum fraction of the parton entering the hard scattering is then 

— = xx TT Zi , (B.3) 
x - LJ - 

ieCi 

where Z; L are the collinear splitting fractions for radiation of secondary partons ki in 
the region C\ collinear to p\. The phase space d& n is then given by 



d® n =dp 2x dp 3x ( Kout - \p2x\ ~ \P3x\ -^\kix\ \ 5 I x B -xx ± Yl 

\ i J \ ieCi 

6 [Plx + Yl kix + kix ) S ( P3x + ^2 k ™ + 7^2 kix j ' 

V ieU i J V ieD i J 



(B.4) 
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where the last two delta functions fix the event plane, as shown in (3.11). 
So the radiation factor X in (3.12) takes the form 

f 1 1 n f d 3 k- f 
Zt,6J = x dx{P T j (Xufl) V T | / 1 / d$ n ■ S s , n . (B.5) 

Jo v n! r=t ^ ^ 7 

To resum the secondary parton emissions we use the factorization structure of the 
soft emission factor Ss, n - The phase space d& n can be factorized by Mellin and 
Fourier transforms. We get 



dN ( xx\ \ f dv e uK ° ut 



[d$ = — f^L(^i) 

J n x B J 27ri\x B J 
x -q r ^MP^ e -H\P«*\-iM Y[ € ( Zi ) U(h) , 



(B.6) 



where 



e( z ) = z N ~ 1 , for ked, 
e(z) = 1 , for k $ Ci, 



(B.7) 



and 



U(k) = u 2 {k x ) = e -"(lfe|-iftfc-) f or ky > , rj k < Vo , 
U(k) = u 3 (k x ) = e -"(lfel-iA*x) for ^ < o , ^ < Vo , (B.8) 
C/(A;) = u (k x ) ee e *f 0fe+^3)*« for ^ > ^ . 

The last source wo(fc) corresponds to partons emitted in the beam region, which 
contributes only to the momentum conservation and not to the observable K ont . 
The near-to-planar region K out <C Q corresponds to the region vQ 3> 1 in the Mellin 
variable. 

Resummation can be now performed and one obtains 



n- ,t, <■ ^ \ fdue uK ont f 1 dx! f dN fxxA N _ , 



11 J 2tt 

a=2,3 J 

where £5 is the full radiator for the secondary emission and is given by 

/d 3 k 
—W s (k) [l-U{k)e{z)\ . (B.IO) 

The 1 in the brackets here represents virtual gluon emission. The radiator £g depends 
on the source variables, v, f3 a , N, rjo, on the primary parton momenta p a , which are 
given in terms of Q,x,£ and the recoil components p ax (the other components can 



vdPadpax -v(\ Pax \-if3 aPax ) -£ s 



(B.9) 
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be neglected in the soft limit), and on the subtraction scale \x needed to regularize 
the collinear divergence. Here Ws(k) is the two- loop distribution for the emission of 
the soft gluon from the three hard partons p a in the configuration 5: 



Ws(h) = y (wsa(k) + w Sb (k) - -^w ab (k)^j , a ^ b ^ 5 , (B.ll) 

where w a b(k) is the standard distribution for emission of a soft gluon from the ab- 
dipole 

Wab{k) ~ ^57 ' aM " to) ' { } 

and a s is in the physical scheme [19]. 

B.2 Factorization of incoming parton distribution 

In the present formulation, the factorization (3.13) of the distribution J results by 
splitting the source 

[1 - U{k)e{z)\ = [1 - U{k)\ + [1 - e{z)\ U{k) , (B.13) 

so that the radiator £ can be expressed as the sum of two terms (we write explicitly 
only the N and \i dependence) 

S s (N, f i)=n s + r s (N,»), (B.14) 

with 

d 3 k 



— W s (k) [1 - U(k)\ , 

TTLd 

r s (N,n)= [ —Ws{k) [\-z N - 1 ] U(k). 

JCi km 



(B.15) 



f d 3 k 

'C-i 

These two terms will be evaluated for small K out ~ z/ _1 . 

Consider first the collinear singular piece Tg(N, fi). Here the source U (k) provides 
an upper bound for the integration frequencies of order v^ 1 ~ K out . This is shown 
by using, in the integral (B.15), 



[1 - u a {k x )\ ~ 6 ( \k x \ - p- 1 ) , Pa = e^u^l + M , (B.16) 

for a = 2, 3, which is valid for large v within SL accuracy (see for instance Appendix C 
of Ref. [11]). Using this approximate expression one also shows [12] that Ts(N,u) 
does not depend on i]q, as long as K out is taken in the range (2.11). As shown in 
[12], r(iV, v) is the soft part of the anomalous dimension that evolves the parton 
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distributions V from the small subtraction scale to the scale p^ 1 . It is accurate 
at two loop order provided we use the coupling in the physical scheme [19], and 
it is diagonal in the configuration index since soft radiation is universal and does 
not change the nature of the incoming parton. But in general one needs the full 
anomalous dimension (including the hard non diagonal pieces) so that the incoming 
partons may change from a quark to a gluon and vice versa. The resulting exponent 
T becomes a matrix. Upon integration over x\ and the Mellin variable N 

r'dxi rdN^y e _ T(M K „(*b \ (B17) 

Jo xi J 2m \ x B J V x / 

one obtains the parton distribution V(K ont ) evolved from the subtraction to the hard 
scale K ovA , up to terms that are beyond SL 1 . The \x dependence in V(/i) and T(N, /i) 
is cancelled. We do not consider here power corrections since they are of second 
order (1/Q 2 ), see [4]. 

We consider then the piece TZs- This is a CIS quantity ([1— U(k)] — >0 for k t — >0) 
which produces the radiation factor As in (3.13) 

A s (K out , x, £, Q, r] ) = / As{v, x, £, Q, r] ) , (B.18) 

with the Mellin moment 

As= I n ud/3adpa * e -»(\p a *\-if>a Pa *) e -Ks ( B19 ) 

J a=2,3 ^ 

All the PT contributions are finite and they can be evaluated to SL accuracy by 
using the approximation (B.16) for the sources. The PT resummed expression 7Z^ T 
is then obtained by evaluating the integral (B.15) of TZs using (B.16). 

However the virtual momentum kt, entering the SL resummed running coupling 
in TZs, cannot be prevented from going to zero. Although the contribution from this 
NP region is (formally) highly sub leading (power suppressed), it is phenomenologi- 
cally quite important for K out not too large [3]. For k t very small (in particular for 
\k x \ v~ x ~ -Kout) the approximation (B.16) is no longer valid and we have instead 

[l-u a (h)]~v(\k x \-ip a k x ), (B.20) 

so that the leading NP correction to TZs is 

5Ks~u- 5K { 5 J t , (B.21) 

corresponding to a shift in K out in (B.18). To evaluate this NP shift we use the 
dispersive method [4] to represent the running coupling and we evaluate the leading 
power-suppressed contribution from (B.20). 



L The more accurate expression \jv{p 2 1 )'P(P3 1 ) is required for some of the NP contributions. 
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The two contributions to the radiator 



K s = 7^ T + SKs , (B.22) 

together with the expression for the radiation factor As will be described in the next 
two subsections. 

B.3 PT Radiation Factor 

The PT radiator TZ$ T to SL accuracy coincides with a particular case of the PT 
radiator computed for T m in e + e~ annihilation (where we set 7 = f3\ = in eq. (4.8) 
of Ref. [10]). For K out in the range (2.11) this contribution is r] independent, see [12]. 
The explicit calculation is given in Appendix C, and the result may be conveniently 
separated into two terms: 

nr^R ( S(p2) + R i S( P s): (B.23) 

where 

RfliPa) = \C?r{ Pa , C?QV) + Cf V(Pa, CW) , (B.24) 

with the function r given in (3.22) and colour charges Cjjp and hard scales given 

in (3.21). The term corresponds to the emission in the up-region off parton p\ 

>(*) 
43 



and p 2 . Similarly for For the PT evaluation we can expand to SL accuracy 



KF = Rs (u, x, £ Q) + [c TlE + C<? In sjl+il + Cf In y/l+^j r'(u, Q) , 

(B.25) 

with Rs the DL exponent function introduced in (3.20) and r' the SL function (3.23). 
Since the radiator is independent of the recoils we can integrate over p ax and, using 
the expansion (B.25), we obtain the Mellin moment 



o(<5) 
- e - R la 



r y (Cg> r') T (C$ r>) , r' = r>, Q) , 



^ (B.26) 



~ e - R s( u > x £>Q) e -lEC T r' 



with T given in (3.25). 

The PT radiation factor is obtained by integrating over the Mellin variable v. 
We make use of the operator identity 



/ 



dv e vK ovA, I 

Giy) = — t — -v G(K£) (B.27) 



2tiiu r [ 1 + a 
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for any logarithmically varying function G. (To prove this, multiply both sides by 
the T-function operator and use the definition T(z) = J °° dx x z ~ 1 e~ x .) Thus to SL 
accuracy we have 



^ T (^o U t, x, Q) = ^^A^ , (B-28) 

T {l + C T r'{K ovlt ,Q)) 

which corresponds to the PT result in (3.19), (3.22) and (3.24). 
B.4 Radiation Factor including NP corrections 

The NP correction to the radiator is computed in Appendix D following the stan- 
dard procedure (see [10]): we assume that the running coupling can be defined 
even at small momenta via a dispersion relation [4], and we take into account the 
non-inclusiveness of the observable via the Milan factor. The NP correction is pro- 
portional to the parameter A NP , given in (D.6), expressed in terms of the integral 
of the running coupling in the small momentum region. This parameter takes into 
account merging of NP and PT corrections at two loops and is the same as entering 
all the jet-shape distributions so far studied. We obtain 

SRs = v\ NP (cf ] In + °a ] ln T^t) ' C = 2e" 2 , (B.29) 

V ^ a=2,3 \ Vax \ J 

where are the scales introduced in (3.21). The NP radiator is made up of three 
contributions, one for each emitting parton. The hard scales are determined by large 
angle NP gluons; they are therefore the same as in the PT case. What makes the 
difference between the three contributions is the NP gluon rapidity cutoff. Actually, 
when a NP gluon is emitted from the incoming parton p 1 , its rapidity is bounded 
by the experimental resolution rjo. This is not the case for the emission from p2 or 
P3, where it is the out-of-plane momentum of the recoiling parton which provides an 
effective rapidity cutoff. 

The first order NP correction to the Mellin moment .4,5 in (B.19) is given by 

A s (u) = TT f ud & d P*° e -*(| PM H&* M ) e -< T { i _ 51Zs} > (b.30) 

where R^ T is the PT radiator in (B.23). The first term yields the PT contribution 
in (3.19) and so we can write 

Mv)=Ar(v)-v\ NP fs(v). (B.31) 
After the p ax integration the NP correction is given by 

(B.32) 
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with x(P) = l n yl + P 2 + /3tan _1 /3. This function is growing as ir\/3\/2 at large (3. 
As a consequence, the /3 a integration in (B.32) is no longer fastly convergent and we 
cannot expand the PT radiator as in (B.25) but we need to keep its exact expression. 
The region of large |/3 a |, which corresponds to \p ax \ <C K ontl gives the leading NP 
correction. 

The radiation factor A is obtained by performing the v integration so that 

A s (K out ) = A p s T (K out ) + 8A s (K out ) . (B.33) 

In Appendix E, see also [7, 10], we show that 8 As can be expressed (to first order) 
as a shift (see (3.29)) 

/J, . pvK OVl i 
f 6 (u) ~ ■ d Kont Ar(K ont ) , (B.34) 



where the shift is given by 

A ^ a=2,3 



n (5,NP) r r (5) 
2 A out G T 



(B.35) 



Here the NP hard scale is defined by 

n (5,NP) tn {8) ( M S K>\ /l,M s ) 



~ In + O (r') , K out = e'^K out , 

A out 



(B.36) 



where r' = r'(K ^ t ,Q) and the choice of factors of 2 and e 7B is so as to give the 
simple form for equation (B.42) below. The function is given by 

Ei'HK^) = r ^-<W<><^> ( nf>) 1 , (B.37) 

where the appearence of the parton distributions comes from the mismatch between 
the correct scales p" 1 and the SL approximation K ont in (B.17). They introduce a 
dependence on r, / and xb which we suppress. 

As explained in detail in [10], these functions result from an interplay between 
PT and NP contributions. They can be expressed as the average of the rapidity 
length over the Sudakov factor for emitting in the up- or down-region for a = 2 or 
a = 3 respectively. Finally, 

Cf< = (a? (A- 1 ,) - ) ( - 1 | (B.38) 




29 



is a regular function of r'. 

The shift can be expressed as a sum of partial shifts as in (3.30) with the rapidity 
integrals given by 



,(<*) 



Y (S) _ y(S) _ ^"Qi 
I 2l — 2 31 — 111 q 

^C^, 4" = ^, (B ' 39) 

r 23 ~~ -^3 l-^outj + L-3 , ^32 — -^2 l-^outj T ^2 ? 

where we have introduced the NP functions L a defined as 

4 5) (^out) = In + ^) . (B.40) 

^ -ft out 

These functions have the following asymptotic expansions: from (F.4) one has 

a - ln2 ^ >>1 ' (B41) 

while from (F.5) one obtains the behaviour of L a in the region a s In 2 Q/K out <C 1: 



1 / n( 5 ' PT ) / n( <5 ' PT ) 



(B.42) 



ia s d\nK out QJ 
where the PT hard scale is Qi 5 ' PT) = (^Q^. 

C. The PT radiator 

The PT radiator is given, to SL accuracy, in terms of afr-dipole radiators 

ratfrfafo) = r a >,/3 2 ) + r£(i/,ft) , (C.l) 
where the 'up' and 'down' hemisphere components are given by 

.r,\. . I d3jc v . fu\ ri „-^(lfel-i&fc)l ... ^^_ as ^»M' 

7TC<J 
d 3 k 



r&foft) = / —w ab (k) [1 - e -"(l*.l-<ft*.)] , Wafe (A;) _ , , 
rS(^,A)= / —w ab (k) [1 - e -(IM-^)] , 



(C.2) 
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and k a b,t is the invariant transverse momentum of k with respect to the P a , Pb hard 
partons in (A.l). For the configuration 6 = 1, for instance, we have 

Rr^^a) = y (W, AO +r 13 (u,(3 a ) -± r23 (v,j3 a )) . (C.3) 

To evaluate the afr-dipole radiator r a b(v, f3 a ) we work in the centre of mass system 
of the afr-dipole. We neglect the rapidity cut (2.3): the correction is beyond our 
accuracy [12]. Denoting by P*, P b * and k* the momenta in this system, we introduce 
the Sudakov decomposition 

p; = ^(i,o,o,i), p; = ^(i,o, o,-i), k* = a p* a + pp; + « , (c.4) 

where the scales Q 2 b = 2(P a P b ) are given in (3.3). Here the two-dimensional vector 
k is the transverse momentum orthogonal to the afr-dipole momenta (k 2 = &^ bt ). We 
have then 

a s (n 2 ) d 3 k d 2 nda k 2 „ s 

w a b{k) = — — ^ > = , «>7^- C.5) 

TTK 1TLU IT CX Qa b 

Since, neglecting the recoils, the outgoing momenta P2 and P 3 are in the yz- 
plane, the Lorentz transformation is in the yz-plane and our observable k x remains 
unchanged. However, the boundary between the 'up' and 'down' hemispheres be- 
comes non-trivial, and different for each dipole, so we must take each in turn. 

The 12-radiator with 

k* = aP; + (3P* + k , (C.6) 
has up and down hemispheres given by 

tt ^ QTm n ^ QTm (n _v 

U : K y > —a , D : K y < —a , (C.7) 

with Tm related to £ by (A. 2). The up component is then 



>, &) = f Ql2 C ^ [1 - e-^~^} e( a+ pL) 7 (Ci 

Jo TTK TT J K 2 /Q 2 2 Oi \ QT M ) 



Integrating over a and K y we have 



J-q 12 \kx\ tt \ ^|Kx| 2 2]/e a; | y 



-Q12 ri 

To show this we introduced K y = £ • \k x \ and used 

^ f\i 1 A — rTo^ 1 /nil i\, 



(C.9) 



/ < ^ ( \kx\VTTp) In ( dfi == ~ ^(2(^1) In ^ . (CIO) 

/^^(l+t 2 ) S V |X| y |/y v / TT^ 2|^| 1 ; 



31 



We extended the t- integration to infinity since it is convergent, then we integrated 
over t by expanding a s to second order. Corrections are beyond SL accuracy. 
Finally using (B.16), we obtain 

M , a ,_ n f Ql2 dk x a s (2k x ) ( Q 12 1 QT M 



Similarly the 'down' hemisphere contribution is 

d / ,, > „ [ Q " dK cjlK) f, Qn 1, QT„\ 
^ &) = 2 1- — I" 1 2fc - 2 It J • <C12) 



The 13 radiator has analogous results 



ei a -^- 2a - ( /^L + K - 



and similarly 



(C.13) 



d / ^ f Ql3 dk x a s (2k x ) ( Qis.l, QT M 

The 23-radiator has a slightly different structure: with 

fc* = aP 2 * + /3P 3 * + k , (C.14) 

the up and down hemispheres are given by 

TT QTm (3 — a QT M f3-a 

U:K y > 2 i- 2 r ^ ^ 2 i- 2 r (c - 15) 

The up component is then 

r&frfc) = [ Qh ^-^ f 1 ^ [1 - e -d-l-^-)l 



(C.16) 



QTm y Q 2 T M Q23 
Performing the integrals as above we obtain 

r 23 (v,[3 2 ) = 2 I -—In—, (C.17) 



d/ /o \ o f Q2B dk x a s (2k x ) Q 23 , . 
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Assembling the various dipole contributions and including hard collinear split- 
tings then yields, to SL accuracy, 

1 0=1 (C.19) 
+ 5 (Cf - Cf ) [r(p 2 , QT M ) - r(p 3 , QT M )] . 

Finally, we note that the terms on the second line contribute only at SL level, since 
the DL terms cancel. We may therefore change at will the hard scales in these terms, 
yielding the result (B.23). 



D. NP corrections to the radiator 



We consider the NP correction Sr ab to the afr-dipole radiator. In this case, as we shall 
see, we need to retain both the recoil momenta p2 X and p3 X and the rapidity cut i]q . 
We write the integral in the afr-dipole centre of mass variables a, (3 and k introduced 
in (C.4) and, to obtain the NP correction 5r ab , we perform the following standard 
operations: 

• the running coupling, reconstructed by two loop emission, is represented by 
the dispersive form [4]. Then, the afr-dipole radiation w a b(k) is written in the 
afr-centre of mass system (see (C.5)) in the form 

a B (/c) [°° dm 2 a cfi (m) 

Wab{k) = ^ = l .r^+^r (m) 

• to take into account the emission of soft partons at two loop order [6], we need 
to extend the source u{k x ) to include the mass m of the soft system. We assume 
k x = k cos — > \J k 2 + m 2 cos 0, with the azimuthal angle of k. Similarly we 
introduce the mass in the kinematical relations such as af3 = (n 2 + m 2 )/Q 2 ab 
for the afr-dipole variables; 

• we take the NP part 5a c s(m) of the effective coupling. Since it has support 
only for small m, we take the leading part of the integrand for small k, and m. 
In particular we linearize the source U (k) 



[1 - U(k)] -> vVk 2 + m 2 I cos 0| &(r] - rj k ) . (D.2) 

Recall that r] k is the rapidity of k in the Breit frame (2.1). Here we have 
neglected terms proportional to j3 2 and /3 3 since they vanish, by symmetry, 
upon the f3 a integration; 
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the recoil component p ax of the outgoing parton p a does provide an effective cut 
in the soft gluon rapidity along the outgoing parton [10, 7]. This is due to a real- 
virtual cancellation which takes place when the angle of the outgoing parton 
p a with the event plane exceeds the corresponding angle of the soft gluon. The 
detailed analysis of real and virtual pieces entails that the contribution from 
the observable \/ k 2 + m 2 1 cos0| in the linear expansion of the source (see (D.2)) 
has to be replaced by 

-a\p ax \, (D.3) 



m 2 I cosi 



V k, 2 + m 2 cos + apax 



\" P = 2Mc Kout dm-^pi, c Kmt = ±. (D.5) 



with a the Sudakov variable in the a6-dipole centre of mass in the forward 
region a > ft; 

• to take fully into account effects of non-inclusiveness of jet observables at two- 
loop order, we multiply the radiator by the Milan factor [6, 16] 

3 (1287T + 1287rln2 - 357r 2 )C j4 -57r 2 n f _ . 

= m rn^Tf • (D - 4) 

using rif — 3; 

• the NP correction is finally expressed in terms of the parameter 

5a cS (m) _ 2 

IT OUt 7T 

After merging PT and NP contributions to the observable in a renormalon free 
manner, one has that the distribution is independent of /ii and one obtains 

A NP = c Kout M^j jaoOi/) - a. - (3 ^ (in £ + |- + l) } , (D.6) 
where 

r, „ f Q7 k 2 \ 5 a 11N C 2n f 

a s = a m (Q) , K = C A ^—-—J- -n f , A, = - -f- . (D.7) 

The X factor accounts for the mismatch between the ms and the physical 
scheme [19] and c*o (/•*/) is the integral of the running coupling over the infra- 
red region, see (3.31). 

The numerical coefficient cx out depends on our observable K out . For instance, 
the shift for the r = 1 — T distribution is 

-r-tr) = -7 -(r-A r ) , A T = C F ^ , c T = 2, (D.8) 

dr dr c Kout 

where Cp enters due to the fact the 2-jet system is made of a quark-antiquark 
pair. 
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We recall that these prescriptions correspond to taking into account NP corrections 
at two-loop order in the reconstruction of the (dispersive) running coupling and in the 
non-inclusive nature of the observable. We implement the rapidity cut by expressing 
rjk the soft gluon rapidity in the invariant form (2.2). 

D.l Dipoles 12 and 13 

Consider first the contribution of the la-dipole (a = 2,3). We decompose the gluon 
momenta along P 1 and P a taken in their centre of mass system: 

P; = ^(1,0,0,1), P* = ^(1,0, 0,-1), k* = aP: + (3P* + K. (D.9) 

and we use the expression in (D.3) in our linearized source, so that the recoil momen- 
tum of parton p a provides an effective rapidity cutoff for a gluon emitted collinear to 
P a . In the region in which k* is emitted close to Pi, i.e. (3 ^> a, we can write gluon 
rapidity in terms of the Sudakov variables (D.9) as follows: 



1, (x B P + q)k , Q^/k 2 + m 2 
rih — - m ; — ~ In k , 

lk 2 x B (Pk) aQ\ a x ' 

so that the NP correction to the la-dipole radiator is given by 



Sr la = 



uM 



7Y 



dm 2 Sa e s(m) 



-d 



dn 2 



dm 2 J k 2 + m 2 



ha , 



(D.10) 



ha = J J ^ ( y/iJ+m 2 cos (p+ap ax -a.\p ax \) © (vo + ln ^7 



aQ\ a x 



k 2 + m 2 

(D.ll) 

If we consider the two rapidity cutoffs discussed above we have that the a integration 
is restricted to the region 



a m < a < (Xm , 



\/ k 2 + m? Q e~ 



a M = 



\f k 2 + m 2 



\Pa 



(D.12) 



giving 



ha = J ^~ J ~ ( V^+m 2 cos <fr+ap ax -a\p ax ^j 

2 f~2~\ 2 ( , 1 Q laX i l Qla(\ , 

= — v k + m \ Vo + hi 1" m i r > C — 2e 

V Q \Pax\ J 

We thus obtain the NP correction to the la-dipole radiator 



Sr 



la 



v\ NP I Vo + In 



Q 



+ In 



Pa 



(D.13) 



(D.14) 
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D.2 Dipole 23 

Again we decompose the emitted gluon momentum along P 2 and P 3 

P* = ^(1,0,0,1), P* = ^(1,0, 0,-1), fc* = aP* + /3P 3 * + «. (D.15) 

In this case, if A;* is emitted close to P 2 * (a > (3), its rapidity is cut by the recoil 
component p 2x , while it is p$ x which provides an effective rapidity cutoff when k* is 
close to P 3 * (a < (3). Thus it is convenient to split the radiator into forward (a > (3) 
and backward (a < (3) regions. In the backward region one may relabel f3 as a. 
Performing the substitution in (D.3) in the linearized source we obtain 



Sr 



23 



ja 
1 23 



vM 



7Y 



dm 2 5a e g(m) 



-d 



dK 2 



dm 2 J k 2 + m 2 



(-^23 + ^23) ) 



J 7^ J ^ ( VK 2 +m 2 cos (j)+ap ax -a\p ax \^ Q - 



\J k 2 + m 2 



Q 



23 



(for a = 2, 3), so that the integration limits for a become 

a m < a < a M , ol Ti _ 
giving the following result 

" aM da 



Q 23 



\/ k 2 + m 2 

\Pax\ 



(D.16) 



(D.17) 



23 = / 7T [ — ( V/t 2 +m 2 cos (f)+ap ax -a\p ax \) 

Q23C 



= —V k 2 + m 2 ( In 



7T 



\Pa 



(D.18) 



C = 2e" 



Putting together the two pieces one obtains the NP correction to r 2 3 



Sr 



23 



NP f j n Q23C , ,„ Q23C 



|j> 



+ In 



2.r 



3x I 



(D.19) 



In conclusion, assembling the contributions from the three dipoles, one is left with 
the expression in (B.29). 



E. Evaluating the NP shift 

Here we evaluate the NP shifts SK^ t introduced in (B.34). Using the operator 
identity (B.27) we may write 

__\NP 9 



SA s (K out ) = — , fs(K^) . (E.l) 

-"■out r ( 1 + 3 



dlnK out 
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Applying this differential operator yields 



SA s {K out ) 



-X 



NP 



K out T(l + R') 



n 



a=2,3 J -°° 



d(3 a e 



T(l + PI) 



+ E c « ] {R'HCKutQ^) + + K) + R'x(Pa) - 1} 

a=2,3 



(E.2) 



with X out = K ont e~^ E , E! = C T r'(K~^ t , Q) and R$ defined in (B.24). The final term 
may now be integrated by parts to give 



5A s (K out ) = 



-A 



NP 



R! 



n 



(S) 



K out Til + R')^^ tt(1 + #) 

+ E ^ I MC^Q?) + iKl + #) + \ + + f ) ft tan- 1 (3 a 



a=2,3 



(E.3) 



For all but the final term in the brackets, we may expand the radiator as in 
the PT calculation and perform the integrals over f3 a . The final term is only slowly 
convergent and must be treated with extra care. Thus 



out J 



-5K, 



(S) 



R' 



out ,- A s (K out ) 



K, 



out 



(E.4) 



with 



SK 



(S) 
out 



a=2,3 



111 



CQ'a 



(5) 



K, 



+ ^(1 + C T r') + \i) 1 + 



out 



1 + Cff r' 



+ ^ ) fc 1 t) 



(E.5) 



The term called m '(K-^ t ) arises from this final term. For the term proportional to 
we may expand the radiator in /3 3 and integrate over it, and vice versa, giving 



vr(l + /3 a 2 ) 



A, tan" 1 & . (E.6) 



The reason we may not simply expand the remaining radiator and integrate over 
(5 a is that doing so gives an unphysical result that diverges in the limit r'(K~^ t ) — > 0. 
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This unphysical divergence is regulated by the second derivative of the radiator. We 
write 

/^tan- 1 /? = ^| - ,3 tan" 1 i. (E.7) 

By making the substitution p = K~^ t y/l + (3%: the contribution from the first term 
is written in terms of the function E a , discussed in the next appendix. The second 
contribution gives a fastly convergent integral, so we can expand the radiator and 
integrate as before. We find 

«P<*i> - (i - f ) (*><*3> + ^§^) ■ (-) 

where the E a functions are discussed below. Thus we recover the result (B.35). 



F. The E a functions 

For the functions Ea\K~* t ) we obtain the following result: 

(F.l) 

where R' la , R'[ a and R'" a are the first three logarithmic derivatives of Rf}(p) evaluated 
at p = K~^ t , and the functions N and X are defined by 

N (t) = 4= r dx e ~ x2+t2 > = -^^P- = r dzz e- z ~ 2t ^ . (F.2) 

V 71 " Jt 8 dt Jo 

They have the following asymptotic behaviours 

t>1: JV(f) = ^( 1 -^ + ? + "')' x w = + ° («"*)• 

f«l: JV(t) = l--?t + t 2 -i^L + ..., A'(t) = 1 - + 4i 2 + O (« 3 ) . 

V VI" OV7T 2 

(F.3) 

In the region -R' la / \plW[ a 1, using the expansions in (F.3), we get 



(FA) 
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On the contrary, for R' la / ^2R'{ a <C 1, we obtain, up to contributions O (, 

R' la 1 d\nV(K out ) RZ 



0. 



TT 



2R» a R'{ a 2R'[ a d\nK ont 
Starting from the definition (B.24) we find 



3^ 



+ 



< = ^ s (2X out )Ucfhn^_ 

TT \ 2 K out 



C a ln 



*2. = ;*« S (2^t) + §« s 2 (2K out ) ±cf In 

7T 7T Z \ 2 A ou t 



+ C hi 



2K 



out 



, 2/5 (5) 2 



and therefore 

^(^ut) 



TT 



AS) n {S) AS) n (S) 



TT <91nP(^ ut) , A) 



2K, 



out 



4a s 9 In if, 



+ 



out 



(F.5) 



(F.6) 



(F.7) 
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